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ABSTRACT 

In  a  laboratory  experiment,  turbulent  mixed  regions  were 
generated  in  a  linearly  density-stratified  fluid  and  their  be¬ 
havior  was  studied.  Such  regions  may  occur  in  nature  in  the 
atmosphere  and  in  the  ocean.  Particularly  during  their  early 
,  history,  the  shape  of  such  regions  is  influenced  by  the  inter¬ 
acting  effects  of  turbulence  and  buoyancy,  culminating  in  the 
occurrence  of  a  maximum  thickness  and  subsequent  vertical  col¬ 
lapse.  A  Richardson  number  (equivalent  to  the  ratio  of  the 
characteristic  turbulence  time  and  the  Vaisala  period)  was  found 
satisfactorily  to  correlate  the  data  obtained,  together  with 
those  previously  obtained  by  other  investigators  with  self-pro¬ 
pelled  bodies.  An  estimate  is  made  of  the  degree  of  mixing  that 
takes  place  inside  a  turbulent  mixed  region  during  its  growth  in 
stably-stratified  surroundings;  the  effectiveness  of  this  mixing 
determines  the  ultimate  thickness  to  which  the  mixing  region 
collapses. 

INTRODUCTION 

Stable  density  stratifications  commonly  exist  in  the  at¬ 
mosphere  as  well  as  in  the  ocean,  and  mixed  turbulent  regions 
are  created  in  these  media  in  a  number  of  ways,  for  example, 
mixed  patches  produced  by  over-turning  internal  waves  or  the 
remnants  of  chimney  trails  or  thermals  in  the  atmosphere.  Also 
bodies  traveling  through  such  media  produce  a  quite  different 
turbulent  wake  behind  them  than  in  a  homogeneous  medium.  The 
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present  study  is  concerned  with  the  time  history  of  a  slice  of 
such  a  wake  from  its  generation  onward,  which  is  closely  equi¬ 
valent  to  studying  the  change  of  wake  shape  with  distance  behind 
the  body. 

In  homogeneous  surroundings,  a  turbulent  wake  grows  in  size 
due  to  turbulent  diffusion  at  the  boundary,  the  rate  of  spreading 
being  directly  proportional  to  the  turbulent  intensity.  Due  to 
spreading  and  turbulence  decay,  the  turbulent  intensity  and 
consequently  the  spreading  rate  continuously  decrease  with  time. 

In  a  stratified  fluid,  turbulent  diffusion  produces  a  density 
distribution  inside  the  turbulent  wake  different  from  that  of 
the  surrounding  fluid.  At  the  same  time  buoyancy  forces  inter¬ 
act  with  the  turbulence  causing  changes  in  the  latter,  see  Webster 
(1964),  for  example.  As  the  turbulent  intensity  decreases  and 
the  stabilizing  buoyancy  forces  increase  with  time,  a  condition 
is  reached  wherein  turbulent  diffusion  in  the  vertical  direction 
is  completely  luppressed.  At  this  instant,  the  wake  or  mixed 
region  has  reached  a  maximum  vertical  thickness.  If  the  fluid 
inside  the  region  is  completely  mixed  and  non-turbulent,  all  the 
interior  fluid  will  flow  to  the  level  of  the  center  of  the  wake, 
and  with  time  the  thickness  becomes  thinner  and  thinner,  as  was 
shown  by  Wu  (1969).  With  incomplete  turbulent  mixing,  the  fluid 
inside  will  not  have  a  uniform  density  and  only  a  partial  verti¬ 
cal  collapse  occurs  as  was  previously  observed  by  Schooley  and 
Stewart  (1963)  and  by  Kennedy,  Stockhausen  and  Clark  (1966). 
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Except  for  the  final  stage  the  history  of  a  turbulent  wake 
In  a  density-stratified  medium  is  shown  herein  to  be  governed  by 
the  Richardson  number,  a  dimensionless  combination  of  the  Initial 
turbulent  Intensity  (actually,  the  Initial  spreading  velocity  of 
the  wake),  the  density  stratification  of  the  medium,  gravity,  and 
the  Initial  wake  radius.  The  interesting  parameters  such  as  maxi¬ 
mum  thickness,  time  necessary  to  reach  this  maximum,  and  final 
thickness,  are  all  seen  to  depend  upon  the  value  of  this  Richard¬ 
son  number. 

The  data  for  the  present  study  are  taken  from  experiments 
carried  out  by  van  de  Watering  (1966).  His  experiments  and  pre¬ 
vious  experimental  results  are  reanalyzed  herein.  The  present 
analysis  not  only  explains,  but  even  enables  a  prediction  of  the 
changes  In  cross-sectional  shape  with  time,  which  are  experienced 
by  a  turbulent  mixed  region  in  any  linearly  density-stratified 
medium. 


INITIAL  CONSIDERATIONS 

The  present  study  is  concerned  with  the  history  of  the  cross- 
sectional  shape  of  a  turbulent  wake  (having  an  Initial  radius  rQ 
and  a  spreading  rate  uQ  In  a  homogeneous  medium  cf  density  pQ  and 
viscosity  u)>  produced  by  a  body  which  can  be  characterized  by  a 
diameter  D  and  which  travels  with  a  velocity  V  through  a  density- 
stratified  medium  of  density  pQ  at  the  height  of  the  vehicle,  of 
viacosity  u*  end  density  gradient  dp/d y,  while  the  gravitational 
acceleration  is  g.  Newton’s  law  for  a  fluid  particle  of  volume 
Q  and  density  pQ,  which  is  A p  heavier  than  its  surrounding  fluid, 
states  that  the  gravitational  force,  acting  on  this  particle, Is: 
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■  (na*«)  x  (acc.)  *  (pQQ)  x  (acc. )  [l] 

Pro  re  [1]  the  acceleration,  denoted  by  g#  ,  la  equal  to 

acc.  -  *#  -  C 

For  a  linearly  denalty-atratlfled  fluid  thla  acceleration  la: 


dt* 


I. 


[2] 


where  y  la  the  vertical  dlatance  from  the  particle  to  lta  equi¬ 
librium  poaltlon. 


In  the  preaent  context  rQ  la  choaen  aa  the  length  acale  In  [2]: 


I. 


-r  £*g 
o  or  • 

"o  ' 


-  v« 


From  [2]  It  can  be  derived  that  the  expreaalon  V 

v  pe  ay 

la  the  frequency  at  which  a  dlaplaced  particle,  if  released,  will 
oscillate  around  lta  equilibrium  position  In  a  linearly  strati¬ 
fied  fluid.  It  la  known  aa  the  Valsala -Brunt  frequency.  We  will 
refer  to  it  aa  V 
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In  stratified  fluid  the  following  parameters  are  of  interest: 

the  initial  and  maximum  thickness,  2y  and  2y _ ,  respectively; 

o  max 

the  initial  wake  vertical  spreading  rate,  vqj  the  time  to  reach 
this  maximum  thickness,  tCQl,  and  the  final  wake  thickness,  2yf. 
These  dependent  variables  are  a  function  of  the  previously  listed 
Independent  variables: 


(yo'Vym.x 


-tcci'yf)  ‘  P  (ro*VD'V'>1’8*'po) 


Using  Buckingham's  ir- theorem  one  obtains 


o  vo  ymax  uo  ‘col  yf  | 

__  n 

fro 

u  r  u  p 
o  o  oro 

VUo  '  ro  #  ro  '  ro  J 

*  u 

lD' 

<| 

■p 

'  v*  1 

r  u 

The  ratio's  and  ~  may  be  determined  either  by  direct  measure 
ment  or  by  model  tests  in  wind  or  water  tunnels.  The  determina¬ 
tion  of  these  ratio's  is  beyond  the  scope  of  the  present  paper, 
and  they  are  assuned  to  be  known  in  the  following.  The  third 
independent  dimensionless  number  is  the  Reynolds  number.  Since 
the  flow  is  always  turbulent  the  Reynolds  number  does  not  play  a 
role.  The  five  dimensionless  dependent  variables  are  thus  seen 
to  be  solely  a  function  of  the  fourth  independent  dimensionless 
number,  generally  known  as  the  densimetric  Froude  number.  Be¬ 
cause  of  the  proportionality  of  velocity  of  advance  of  the 
turbulent  front,  uQ,  and  the  turbulent  intensity, 'y  u'2  ,  we  may 
also  see  this  number  as  a  Richardson  number,  defined  in  this 
particular  application  as: 
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For  the  above  reasons  we  may  thus  expect,  that  the  five  dimension¬ 
less  parameters  of  interest  will  vary  only  if  the  Richardson  number 
is  varied.  The  Richardson  number  has  been  used  before  to  scale 
flows  in  which  turbulence  and  gravity  forces  are  interacting  (see 
e.g.  Townsend,  1958,  or  Defant,  1961 ).  It  may  also  be  thought  of 

as  the  ratio  of  a  characteristic  turbulence  time,  r  /u  to  a 

_  o  o 

characteristic  stratification  time,  1/  y  ag. 

Through  multiplication  of  two  dimensionless  numbers  one  ob¬ 
tains  tcolV ag,  which  is  to  be  preferred  above  ^C0iu0/r0  to  n°n- 
dimensionalize  the  time  it  takes  to  reach  the  maximum  vertical 
thickness. 

Next  let  us  see  how  the  laboratory  experiments  prove  that 
indeed  the  pertinent  parameters  associated  with  turbulent  wake 
growth  and  collapse  only  depend  upon  the  value  of  the  Richardson 
number. 


EQUIPMENT  AND  EXPERIMENTAL  TECHNIQUE 

The  spatial  behavior  of  a  turbulent  wake  in  a  density- 
stratified  fluid  is  believed  closely  equivalent  to  the  time 
history  of  a  thin  slice  of  such  a  wake,  considered  as  two-dimen¬ 
sional.  With  this  analogy  in  mind  we  designed  the  experimental 
equipment  and  technique. 


HYDRONAUTICS,  Incorporated 


-7- 

The  experiments  were  carried  out  in  a  transparent  lucite 
tank,  121  x  109  x  31.5  cm,  in  which  a  linear  density  stratifi¬ 
cation  was  obtained  by  mixing  increasing  amounts  of  sodium 
chloride  in  water  and  introducing  these  mixtures  into  the  tank 
in  layers  through  a  diffusor  placed  on  the  bottom  of  the  tank. 

The  layers  were  2.54  cm  thick  and  alternately  colored  with  red 
dye.  In  this  way  24  layers  were  brought  into  the  tank  which 
formed  initially  a  stepped  density  distribution.  Overnight  a 
linear  density  gradiunt  was  obtained  through  molecular  diffusion, 
but  the  dye  remained  in  discrete  layers  because  of  its  relatively 
low  diffusivity. 

At  mid-depth  of  this  linearly  density-stratified  fluid,  a 
spiral  paddle  was  supported  across  the  tank  by  a  0.8  cm  I.D.  brass 
tube  and  rotated  forward  and  back  by  means  of  a  pendulum  type  of 
arrangement  at  the  back  of  the  tank.  Four  rows  of  holes  were 
drilled  along  this  tube.  Just  prior  to  a  test,  a  rod,  heavily 
painted  with  a  blue  dye  which  rapidly  dissolves  in  water,  was 
inserted  in  the  perforated  center  tube.  The  blue  dye  was  used 
to  facilitate  visualization  of  the  turbulent  mixed  region.  The 
above-described  "echnique  to  generate  a  turbulent  mixed  region 
was  first  used  by  Wu  in  1965 .  Two  spiral  paddles  were  used  in 
these  experiments;  one  was  5.08  cm,  the  other  10.16  cm  in  diam¬ 
eter  (see  Figure  l).  With  the  5*08  cm  paddle,  two  different 
pendulums  were  used  so  as  to  vary  the  turbulent  mixing  character¬ 
istics  . 
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The  life  of  the  mixed  region  was  photographed  with  a  16  mm 
movie  camera.  Zero  time  was  taken  at  the  instant  the  paddle  mo¬ 
tion  stopped.  Data  were  taken  from  tracings  of  projected  frames 
of  the  movie  film  (see  Figure  2). 

If  desired, more  information  pertaining  to  the  experimental 
technique  and  apparatus  can  be  obtained  from  the  original  data 
report  by  van  de  Watering  (1966). 

RESULTS 

With  the  three  combinations  of  paddles  and  pendulum,  tests 
in  pure  water  were  carried  out  first.  The  cross-section  of  the 
mixed  region  during  these  tests  remained  more  or  less  circular. 
From  tracings,  the  averaged  (root-mean-square)  radius  was  deter¬ 
mined  and  plotted  versus  time.  The  initial  velocity  of  advance 
of  the  turbulent  front  in  pure  water,  uq,  together  with  the 
initial  averaged  radius  of  the  turbulent  region  rQ,  are  given 
in  Table  1. 

In  density-stratified  surroundings,  the  shape  of  the  mixed 
region  changes  continuously.  Just  after  generation,  the  cross- 
section  is  already  non-circular  and  the  density  inside  the  mixed 
region  is  different  from  its  surrounding  medium,  causing  sta¬ 
bilizing  gravitational  forces  to  counteract  turbulence.  Webster 
(1964)  and  others  clearly  show  that  the  vertical  scale  of 
turbulence  is  reduced.  With  increasing  time  the  cross-section 
becomes  more  ellipse-like  and  eventually  reaches  a  maximum  thick¬ 
ness,  after  which  a  vertical  collapse  takes  place  to  a  finite 
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thlckness,  2yf.  A  series  of  sample  pictures  Is  presented  in 
Figure  2,  In  which  the  above-mentioned  characteristics  are  clearly 
visible.  For  a  typical  time  history  see  also  Figure  3. 

Since  one  of  the  primary  alms  of  the  present  study  was  to 
relate  turbulent  mixing  and  density  gradient  with  the  character¬ 
istics  of  the  mixed  region,  most  attention  was  focused  on  the 
change  with  time  of  the  vertical  ordinate  of  the  boundary  of  the 
mixed  region  at  the  center  of  the  wake,  since  It  Is  this  thick¬ 
ness  which  Is  affected  moat  by  the  interplay  of  turbulence  and 
gravity.  For  several  runs  the  horizontal  extent  of  the  mixed 
region  was  also  analyzed  In  order  to  obtain  information  on  the 
rate  of  collapse.  For  all  tests,  the  actual  measured  values  for 
half  the  maximum  thickness  (y_„),  the  time  It  took  to  reach 

lillX 

this  maximum  (tCQl)  and  half  the  final  thickness  (yf)  are  given 
In  Table  1. 


ANALYSIS  AND  DISCUSSION 

In  a  stably  stratified  fluid,  the  vertical  ordinate  and 
velocity  of  the  boundary  of  the  mixed  region,  both  taken  at  the 
center  of  the  wake  and  given  the  symbols  y  and  v  respectively, 
are  affected  by  the  Interaction  of  turbulence  and  density 
stratification  and  by  the  length  scale.  The  experiments  indi¬ 
cate  that  this  Interaction  is  already  effective  during  the 
generation  stage.  The  data  show,  see  Table  1,  that  at  zero 

time,  y  ,  and  v  are  In  almost  all  tests  smaller  than  their 
o  o 

corresponding  values  in  pure  water,  rQ  and  uq.  In  Figure  4  the 
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d  linens  Ion  less  initial  vertical  ordinate  and  velocity  are  plotted 

versus  the  Richardson  number.  We  have  Included  the  data  points 

from  the  two  studies  with  self-propelled  bodies,  the  experiments 

by  Schooley  and  Stewart  (19^3)  and  by  Kennedy,  Stockhausen  and 

Clark  (1966).  Prom  the  former  study,  the  values  for  rQ  and  uq 

were  obtained  from  the  original  movie  film  of  these  authors, 

l.e.,  from  tracings  of  projected  frames  of  their  film  of  the  two 

runs  In  pure  water.  At  this  point,  It  may  be  mentioned  that  In 

their  analysis,  a  small  error  was  made  In  the  calculations  of 

the  density  gradient.  The  corrected  value  for  -^/ag"  la  2. 19  sec*1 

while  u  was  found  to  be  5.87  cm/sec  and  r  equal  to  O.568  cm. 
o  o 

Por  the  experiment  carried  out  by  Kennedy,  et  al,  the  correct 

Valsale-Brunt  frequency  la:  V  ®6  *  0.25  sec*1.  The  values  for 

u  and  r  were  obtained  from  the  curve  of  the  horizontal  width 
o  o 

of  the  wake,  which  Is  not  significantly  affected  by  gravitational 
effects  until  collapse  starts.  The  values  obtained  from  these 
studies  are  tabulated  In  Table  1. 

The  velocity  seems  much  more  affected  by  the  Interaction 
than  the  thickness,  see  Pigure  4.  The  scatter  of  the  data  does 
not  Justify  a  definite  relation,  but  the  tendencies  are  clear. 

Por  Increasing  Richardson  numbers,  (corresponding  to  a  large 
stratification),  both  y  /r  and  v  /u  tend  to  be  reduced,  while 

OOOO 

for  small  Richardson  numbers  (turbulence  is  large  compared  to 
stratification),  in  the  range  from  0  to  1,  there  Is  no  signifi¬ 
cant  effect.  In  the  ocean,  R^  may  vary  from  about  0.06  to  0.9 
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and  we  may  thus  safely  say  that  in  this  range  the  effect  of 
stratification  upon  initial  conditions  is  negligible  for  all 
practical  applications. 

During  the  growth  stage  the  buoyancy  forces  increase  while 

the  turbulence  decreases  in  intensity.  The  increasing  dominance 

of  buoyancy  culminates  at  some  time  after  generation,  t  in  a 

maximum  wake  thickness  (*2y  ).  The  maximum  thickness,  y  /r  , 

'  max7  max'  o 

where  r  is  the  radius  of  the  circular  mixed  region  in  uniform 
o 

surroundings  at  time  zero,  is  plotted  in  Figure  5a.  In  all 

laboratory  experiments  on  turbulent  wake  collapse,  the  maximum 

thickness  can  be  observed  and  measured  rather  accurately.  Due 

to  the  nature  of  the  variation  of  thickness  with  time  it  is  much 

more  difficult  accurately  to  determine  the  time  at  which  this 

maximum  occurs.  Therefore,  we  may  expect  more  scatter  in  the 

data  of  t  ,  as  compared  to  those  of  y  .  The  data  are  pre- 

sented  in  Figure  5b.  It  is  not  possible  accurately  to  determine 

the  time  at  which  collapse  started  from  the  experimental  results 

as  given  by  Kennedy  and  his  co-workers.  Therefore,  the  possible 

range  of  t  ,  V  ag  is  indicated, 
coi 

During  the  growth  stage,  the  turbulent  mixing  inside  the 
wake  will  produce  therein  a  density  gradient  less  pronounced 
than  that  of  the  surrounding  fluid.  After  the  wake  grows  to  its 
maximum  thickness,  the  fluid  particles  within  seek  their  equi¬ 
librium  level  under  the  action  of  hydrostatic  pressures.  The 
mixed  region  decreases  in  thickness  ar.d  spreads  rapidly  hori¬ 
zontally.  This  phenomenon  is  known  as  wake  collapse  (Schooley 
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and  Stewart,  1963  and  Schooley,  1967).  If  the  mixed  region  is 
non-turbulent  and  fully  mixed,  a  long  thin  wedge  will  be  formed, 
which  becomes  thinner  with  increasing  time;  such  well  mixed 
regions  have  been  observed  and  studied  by  Wu  (1969).  Inside  the 
usual  turbulent  mixed  region  in  a  stratified  fluid,  turbulent 
diffusion  does  not  succeed  in  mixing  completely  the  wake  contents, 
and  it  is  for  this  reason  that  there  exists  an  effective  density 
distribution  within  the  wake  at  the  time  when  the  maximum  thick¬ 
ness  is  reached.  This  density  distribution  inside  the  mixed 
region  should  again  depend  only  upon  the  Richardson  number. 

Wake  collapse  is  driven  by  the  hydrostatic  pressure  result¬ 
ing  from  the  differences  between  the  exterior  (ambient)  and  in¬ 
terior  density  gradients.  The  interior  gradient  is  continually 
increased  as  a  result  of  the  vertical  shrinking  of  the  wake.  The 
mixed  region  thus  reaches  a  final  thickness,  (2yf)),  corresponding 
to  which  these  density  gradients  are  of  equal  value.  The  mea¬ 
sured  values  of  yf  are  given  in  Table  1  and  plotted  dimension- 
lessly  in  Figure  5c, 

The  trend  of  the  data,  shown  in  Figures  5a,  5b  and  5c,  is 

readily  explained.  A  turbulent  mixed  region  in  homogeneous 

fluid  (R^  =  0)  will  continue  to  grow  indefinitely  in  size  and 

y  ,  t  1  and  y-  are  unbounded  at  this  limit.  On  the  other 
max  col  f 

hand  the  non-turbulent  mixed  region  (Wu,  1969)  in  density- 
stratified  fluid  (uq  =  0  and  R^  ♦  «)  is  the  other  limit  for 
which  we  know  that  there  is  no  growth,  i.e.,  t  ^  =  0,  while 
yf  approaches  zero.  In  between  these  limits  the  data  show  a 
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gradual  increase  in  the  value  of  the  parameters  as  the  Richardson 
number  decreases,  i.e.  as  the  wake  becomes  relatively  more  turbu¬ 
lent.  The  opposite  case,  i.e.  decreasing  maximum  thickness,  time 
to  reach  maximum  thickness,  and  final  thickness,  is  observed  when 
the  density  gradient  or  R.^  is  increased.  This  interpretation  also 
explains  why  it  was  observed  that  with  the  larger  paddle  (larger 
R^  the  turbulent  mixed  region  did  not  seem  to  have  much  strength 
i.e.  the  maximum  thickness  was  small  compared  with  the  initial 
diameter  in  a  homogeneous  medium. 


In  the  past,  the  time  t  ^  has  generally  been  considered  as 
solely  dependent  upon  the  density  gradient  (see  for  instance 
Schooley,  1967  and  1968),  but  the  present  analysis  shows  that  it 
is  the  Richardson  number  that  actually  governs  the  value  of 


Let  us  suppose  as  an  approximation  that  at  t 

exists  inside  the  mixed 


t  a 
col 


linear  density  gradient  (dp/dy )lnternal 

region,  see  Figure  6.  The  difference  in  density  of  the  fluid 
at  the  top  of  the  wake,  (l)  in  Figure  6,  compared  to  fluid  in  the 
center  (2),  is: 


P»  -  P°  “  ymax  ’  (dP/3y) internal 

If  we  assume  that  no  further  mixing  takes  place  after  t  ^ 
then  finally: 


yraax  (SP^y)lnternal  '  yf  <3P^y)  external 
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y 

‘'max 


M 

'  'internal  at  t  , 
_ col 

§£■ 

y /external 


[3] 


The  right  hand  side  of  this  expression  is  a  measure  for  the 
degree  of  mixing,  $,  which  takes  place  during  growth  and  which 
we  define  as: 


♦ 


1$ 

v  y' internal  at  t 


1  - 


‘I 


col 


external 


[4] 


expressed  in  percent.  Combining  [3]  and  [4],  the  degree  of 
mixing  is: 


♦ 


y 

max 


x  100 


[5] 


In  Figure  7  the  values  of  yf/ymax  and  ♦  are  plotted  versus  the 
Richardson  number.  Increasing  the  density  gradient  (i.e.,  in¬ 
creasing  the  value  of  the  Richardson  number)  will  result  in  less 
mixing  and  from  Figure  7  it  is  seen  that  indeed  the  degree  of 
mixing  decreases  with  increasing  Richardson  number.  Mixing  is 
generally  poor:  the  degree  of  mixing,  ♦,  is  about  40  percent 
at  the  smallest  Richardson  numbers  and  decreases  if  R1  increases. 
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As  far  as  horizontal  wake  spreading  Is  concerned.  It  should 
be  realized  first  that  a  turbulent  wake  expands  naturally  unless 
opposed,  as  by  hydrostatic  pressures.  While  stratification  acts 
against  expansion  In  the  vertical  direction,  it  actually  enhances 
the  horizontal  expansion.  Initially  the  horizontal  expansion  Is 
dominated  by  turbulence  and  the  effect  of  stratification  Is  small. 
Consequently  the  horizontal  expansion  until  collapse  sets  In,  will 
be  close  to  the  pure  water  case.  However,  once,  the  maximum 
vertical  thickness  is  reached,  the  wake  collapses  vertically  and 
the  fluid  spreads  out  horizontally  under  the  action  of  hydro¬ 
static  pressures.  We  may  thus  not  be  surprised  to  observe  a 
sudden  surge  in  horizontal  expansion  some  short  time  after  col¬ 
lapse  sets  In.  At  this  time  gravitational  forces  start  to  play 
a  completely  dominant  role.  For  several  runs  the  horizontal 
spreading  is  plotted  dimenslonlessly  In  Figure  8.  A  distinct 
discontinuity  in  the  x-t  curves  is  noticeable  and  it  occurs  in 
each  case  a  short  time  after  the  maximum  thickness  is  reached. 

In  experiments  with  a  non-turbulent  mixed  region,  the  area 
of  mixed  fluid  remains  constant  throughout  a  test.  In  a  turbu¬ 
lent  mixed  region,  an  increase  in  cross-sectional  area  with  time 
due  to  turbulent  entrainment  is  observed.  In  stratified  fluid, 
we  observed  turbulent  entrainment  even  after  the  maximum  thick¬ 
ness  was  reached.  A  continuous  Increase  in  area  was  measured, 
although  the  rate  of  area  Increase  decreased  gradually  with  time. 
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The  vertical  collapse  of  the  mixed  region  generated  internal 
waves  in  the  surrounding  fluid  due  to  displacement  of  mixed  fluid 
by  the  ambient  medium.  In  a  turbulent  mixed  region,  even  after 
the  onset  of  collapse,  there  is  still  turbulence,  albeit  that  the 
motion  is  dominated  by  gravitational  effects.  This  turbulence, 
which  tends  to  expand  the  mixed  region,  tends  to  soften  wake  col¬ 
lapse.  The  net  result  is  to  Increase  the  characteristic  time  of 
wake  collapse  over  that  observed  in  the  case  of  the  collapse  of 
a  non-turbulent  mixed  region.  This  time  is  further  Increased 
because  of  the  density  gradient  existing  within  the  wake.  Further¬ 
more,  in  the  turbulent  case  the  mixed  region  typically  collapses 
to  a  final  thickness  which  la  a  large  fraction  of  the  maximum 
thickness.  This  la  shown  in  Figure  7.  The  punch  action  of  the 
impulsive  collapse  la  therefore  not  only  slower  in  time,  but  also 
of  a  smaller  amplitude.  Collapsing  turbulent  mixed  regions  are 
consequently  less  effective  generators  of  internal  waves  than 
non-turbulent  completely  mixed  regions.  All  of  these  effects 
were  observed  in  the  present  experiments. 

Two  other  interesting  features  were  also  observed  which 
seem  to  be  characteristic  for  a  collapsing  turbulent  mixed  region. 
After  reaching  its  maximum  thickness,  mixed  fluid  spreads  out 
horizontally  and  the  upper  and  lower  wake  boundaries  become  more 
or  less  flat;  this  was  particularly  clear  near  and  after  the 
time  when  the  final  thickness  was  reached.  Simultaneously, 

"finger- like”  puffs  of  mixed  fluid  began  protruding  horizontally 
away  from  the  bulk  of  the  wake.  Presumably  these  phenomena  are 
due  to  the  action  in  the  large  turbulent  eddies,  which  are  sup¬ 
pressed  and  flattened  by  the  stabilizing  buoyancy  forces. 
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CONCLUSIONS 

A  turbulent  mixed  region  In  a  density-stratified  fluid 
undergoes  changes  as  a  result  of  the  interplay  between  turbulence 
and  stabilizing  buoyancy  forces.  Some  time  after  generation  the 
mixed  region  reaches  a  maximum  vertical  thickness, after  which 
vertical  collapse  and  horizontal  spreading  take  place.  Due  to 
incomplete  mixing,  the  density  of  the  fluid  inside  the  region 
la  not  uniform,  and,  consequently,  the  mixed  fluid  reache-  an  equi 
librium  shape  with  a  finite  vertical  thickness.  Based  on  the  pre¬ 
sent  experimental  results,  obtained  from  tests  using  paddle  mixers 
as  generators  of  turbulent  mixed  regions,  and  on  other  available 
data  from  previous  studies  by  other  investigators,  the  maximum 
thickness,  the  time  it  takes  to  reach  this  maximum,  and  the  final 
thickness,  are  shown  to  depend  upon  the  ratio  of  the  character¬ 
istic  turbulence  time  to  the  Valsala-Brunt  period. 

This  ratio  is  in  the  form  of  a  Richardson  number.  The 
present  correlation  provides  an  explanation  for  the  higher  values 
of  the  maximum  thickness  and  the  time  to  reach  this  maximum,  as 
obtained  by  other  investigators  with  self-propelled  bodies  under 
conditions  of  relatively  small  Richardson  numbers. 

A  rough  picture  of  turbulent  wake  collapse  is  presented 
with  which,  together  with  a  knowledge  of  the  ratio  of  final  to 
maximum  thickness  (which  can  be  estimated  from  the  present  cor¬ 
relation),  the  degree  of  mixing  inside  the  wake  can  be  estimated. 
Nixing  is  generally  poor,  in  the  order  of  10  to  4o  percent  (in 
a  completely  mixed  fluid  the  mixing  is  100  percent). 


HYDRONAUTICS,  Incorporated 


-18- 

The  collapse  of  a  turbulent  mixed  region  in  a  density- 
stratified  medium  generates  less  pronounced  internal  waves  in 
comparison  with  a  non- turbulent,  completely  mixed  region. 

The  present  experiments  revealed  the  formation  of  a  flat 
top  and  bottom  of  the  mixed  region  and  of  "finger-like"  puffs 
of  mixed  fluid  spreading  horizontally  outward. 
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FIGURE  1  -  SPIRAL  PADDLES;  10.16  cm  AND  5.08  cm  IN  DIAMETER 
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FIGURE  2  -  SAMPLE  PICTURES  OF  TURWLENT  MIXED  REGION  IN  STRATIFIED  FLUID 
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FIGURE  3  -  TYPICAL  HISTORY  OF  A  TUR8ULENT  MIXED  REGION  IN  A 
DENSITY  -  STRATIFIED  MEDIUM. 
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FIGURE  4  -  DEPENDENCE  OF  INITIAL  CONDITIONS  UPON  RICHARDSON  NUMBER. 
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FIGURE  5  -  CORRELATION  OF  MAXIMUM  THICKNESS,  THE  TIME  MAXIMUM 
THICKNESS  IS  REACHED  AND  THE  FINAL  THICKNESS  WITH  THE 
RICHARDSON  NUMIER. 
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FIGURE  6  -  GRAPHICAL  DETERMINATION  OF  FINAL  THICKNESS 
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FIGURE  7  -  CORRELATION  OF  MAXIMUM  THICKNESS,  FINAL*  THICKNESS  AND  DEGREE 
OF  MIXING  WITH  THE  RICHARDSON  NUMIER. 


FIGURE  8  -  HORIZONTAL  SftCAOfNG  REFORE  AND  AFTER  ONSET  OF  COUAFSE 
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